Ring-shaped luminescence patterns in a locally photoexcited electron-hole bilayer 
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We report the results of molecular dynamics simulation of a spatiotemporal evolution of the 
locally photoexcited electrons and holes localized in two separate layers. It is shown that the ring- 
shaped spatial pattern of luminescence forms due to the strong in-layer Coulomb interaction at 
high photoexcitation power. In addition, the results predict (i) stationary spatial oscillations of the 
electron density in quasi one-dimensional case and (ii) dynamical phase transition in the expansion 
of two-dimensional electron cloud when threshold electron concentration is reached. A possible 
reason of the oscillations and a theoretical interpretation of the transition are suggested. 
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Exciton systems in semiconductor nanostructures can 
provide an unique opportunity to obtain nonequilibrium 
superfluid transition in solids [H-Q. Indeed, a remark- 
able phenomenon has been experimentally discovered in 
the system of interwell excitons in GaAs/AlGaAs dou- 
ble quantum wells (QWs) local off-resonant pho- 
toexcitation of electrons (e) and holes (h) gives rise to 
a macroscopic ring-shaped pattern of the luminescence 
spatial distribution when the excitation power is above 
some threshold value. The similar effects were observed 
later in other systems 

The first theoretical explanation of the phenomenon 
was suggested [§- 11 1 within the diffusive transport model 
(DTM). However, the model is not applicable to the ex- 
periments [f| for the derivation of the ring radius [l2j . In 
particular, the ring radius R is derived within the DTM 



R = Aexp (— 2irD e no/ P e: 



(1) 



where D e is the diffusion coefficient for electrons, no is 
the equilibrium two-dimensional (2D) electron density, 
P ex is the stationary photoexcitation power, and the elec- 
tron depletion length A ^> R 9] . According to the model 

0, the diffusion of holes is the only reason why they 
move out of the laser excitation spot. However, accord- 
ing Eq.([T]) the ring radius does not depend on the hole 
diffusion coefficient Dh (in [t| Dh is in the exponent also), 

1. e. the ring forms even if all holes are left into the ex- 
citation spot {Dh = 0). This clearly unphysical result 
is not a consequence of the limiting case A 3> R, but is 
valid in the general case [ll|. Moreover, the DTM can 
not explain in principle the condition on the excitation 
power, i.e. that the external ring appears only when the 
power exceeds some critical value, P ex > (P e x) c ~ 250 
/iW So the development of a different theory of the 
ring pattern formation is necessary. 

Here we present both quasi one-dimensional (ID) and 
2D simulations of the spatiotemporal dynamics of opti- 
cally generated electrons and holes that are located in 
two coupled layers (Fig.l). We show that due to (i) large 



in-layer Coulomb interaction in the pumping region at 
high excitation intensities and (ii) essential difference in 
mobilities and effective masses for electrons and holes in 
GaAs, a macroscopic spatial separation of the electron 
and hole densities appears naturally. In contrast to the 
previous theoretical studies , the in-plane distribu- 

tions for electrons and holes are found self-consistently. 
The distribution of holes is located much closer to the 
excitation spot than the electron one, and their narrow 
overlap gives large ring-shaped luminescence pattern. We 
also find that the dependence of exciton formation rate on 
the electron-hole relative velocity is crucial for the ring- 
shaped pattern formation. More surprisingly, we have 
clearly observed stationary spatial oscillations of the elec- 
tron density outside the luminescence ring in ID case. A 
possible origin of the oscillations is a kinetic instability of 
two interacting electron counter flows: faster flow from 
the center to the periphery and slower back flow from the 
periphery to the luminescence region due to the depletion 
of electron density there. 
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FIG. 1: Schematic of the optically excited electron- 
hole bilayer. Both stationary laser pumping and lumines- 
cence are perpendicular to the layers. The charge separation 
between the layers is due to the external gate voltage V ga te- 
The interlayer distance d is essentially increased for the sake 
of clarity. 

Model. — To describe the dynamics of N hot electrons 
and holes, we use the Newton equations of motion com- 
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bined with the exciton formation condition. The dimen- 
sionless equations read 



r. = 



d 2 



3/2 ' 



(2) 



r h + cir 



1 = E 



r [(r* h - r^)2 + <*> 



3/2 ' 



where vectors r* and are in-plane positions of i-th 
electron and j-th hole (1 < i, j < N); d is the interlayer 
distance. The left-hand side of Eqs.([2]) accounts electron 
and hole inertia and momentum damping (mainly due 
to interaction with acoustic phonons), while the right- 
hand side describes the in-layer electron-electron and 
hole-hole Coulomb repulsion as well as the interlayer 
electron-hole attraction. Hereafter, we normalized time 
by t s = y/emln e /e and all distances by r s = %/m%n1, 
where e and \i e are electron charge and mobility, respec- 
tively, m* is effective electron mass and e is dielectric 
constant of the layers. The coefficients in Eqs.([5]) are 
c\ = m* e ^L e / {m* h ^ih) and C2 = m* e /m* h . To estimate the 
parameters, we use well-known experimental values for 
high-quality undoped GaAs/AlGaAs QWs 13J getting 
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1, c 2 ~ 0.1, t s 



10" 



10 cm that gives 



10 for interlayer spacing of about 10 cm [14 1 



Stationary optical excitation of electrons and holes was 
modeled by generating the particles in random positions 
inside the excitation spot of radius tq with rate p. Initial 
velocities of electrons and holes were also chosen ran- 
domly within intervals (— V e , V e ) and (— Vh,Vh) with ra- 
tio Vh/Ve — C2 to allow the momentum conservation for 
each e-h pair created. During the dynamical evolution of 
the carriers, an exciton formation happened if an electron 
and a hole were close enough to each other, |r e — < a, 
where a(d) is a phcnomenological in-layer exciton radius 
[lit . an d their relative velocity V = |r e — r h \ was smaller 
than some critical value V c Jig]. To simplify the calcu- 
lations, we did not consider the exciton dynamics and 
neglected the interaction of charge carriers with dipole 
excitons. It means that as soon as an electron and a hole 
had formed an exciton, their dynamics was no longer 
considered and the position of the formation event was 
recorded as a position of photon emission. In addition, 
due to the inevitable restrictions in computational power 
we were able to simulate TV < 10 3 interacting particles. 
For this reason we had to increase the value of a. 

Results for ID case. — The results of extensive ID sim- 
ulations of Eqs.Q combined with the exciton formation 
condition are shown in Fig. 2. The stationary lumines- 
cence density n\ um (x) has a large distant peak that corre- 
lates with the ring-shaped pattern observed experimen- 
tally [j| @, 0, [13. It is remarkable that the smaller 
critical relative velocity V c for the exciton formation, the 



larger both the peak amplitude and the distance between 
the peak position and the excitation spot (— xo,Xq) in 
the center (Fig. 2a). Indeed, at extremely large V c = 00, 
when an exciton is formed as soon as the instant dis- 
tance between an electron and a hole is small enough, 
the luminescence decreases monotonically from the exci- 
tation spot. For the finite but still high critical velocities 
a broad plateau-like maximum appears out of the centre 
in the luminescence pattern. Finally, a narrow high peak 
of luminescence located on the external plateau edge far 
away from the excitation spot develops at lower V c . To 
uncover the origin of this peak, we plot density distri- 
butions for electrons and holes (Fig. 2b). It is seen that 
the hole distribution is located much closer to the ex- 
citation spot than the electron one. The luminescence 
peak forms in the region where electron and hole den- 
sities overlap. In addition, the electron density exhibits 
pronounced spatial oscillations beyond the luminescence 
peak (Fig. 2b). They depend on optical pumping param- 
eters: e.g., the characteristic period of the oscillations 
is proportional to p -1 / 4 , where p is pumping rate. By 
turning off interlayer electron-hole interaction, we have 
concluded that the oscillations originate due to the in- 
layer electron-electron interaction. In particular, they 
can be caused by a kinetic instability of two interacting 
electron counter flows: the flow of hot electrons expand- 
ing from the center meets the returning flow of electrons 
from the steady electron cloud at the periphery. The 
back flow can originate from the depletion of the density 
of slow electrons in the luminescence region due to the 
exciton formation. Nevertheless, a detailed understand- 
ing of the underlying physical mechanism of the electron 
density oscillations still requires further studies. The de- 
pendence of the luminescence peak on the electron-hole 
creation rate p (inset in Fig. 2a) has shown that though 
the peak height increases with p, its position is almost in- 
dependent on p for ID case. This behaviour differs from 
that observed experimentally for 2D samples [5j, where 
the luminescence ring radius grows essentially with the 
increase of pumping intensity. The difference poses a 
question of modeling the full- value 2D expansion of pho- 
toexcited electron-hole plasma. 

Results for 2D case. — Due to the quadratic increase of 
the number of particles, in 2D case without a boundary 
(Fig. 3) the stationary state was not reached for the ex- 
pectation time, though a clear tendency to slow down 
was observed for the expansion of the electron distribu- 
tion. Remarkably, however, that the luminescence ring- 
shaped pattern forms dynamically at t/t s > 300, fastly 
expanding afterwards (Fig. 3a) (cp. |17j). Moreover, the 
spatiotemporal dynamics of the 2D electron density (not 
shown) has revealed very nontrivial details: at t/t s ps 300 
the electron cloud ejects a long "tail" towards the outside 
of the excitation area. Due to the corresponding deple- 
tion of the internal part of the cloud, the in-layer sepa- 
ration between electron and hole distributions becomes 
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more pronounced giving rise to a sharp luminescence ring 
in the bilayer plane. 

To roughly understand why the dynamical transition 
occurs in 2D case, let us reduce Eqs.([2]) to the case 
when the carrier distributions can be described only by 
collective ring coordinates r$ = N^ 1 J rrii(r)d 2 r, where 
i — e,h. In the limit r e 3> r/j the equation of motion for 
the electron ring reads 



N N 
— ln(2r e /A)--, 



(3) 



where A <C r e is a width of the electron ring. (For sim- 
plicity, we neglect the dependences N (t) and A(r e ).) 
The first term in the right-hand side describes the inter- 
action of electrons within the ring, while the second term 
originates from the electron-hole interaction. The loga- 
rithmic factor in the electron-electron interaction term 
is absent in ID case. It results in a fast expansion of 
the electron ring for ln(2r e /A) > 7r, i.e., when r e has 
amounted to some critical value. 

We have confirmed that the principle of the ring pat- 
tern formation in 2D case is the same as that in quasi 
ID case: the luminescence peak appears at the narrow 
overlap of separated electron and hole density distribu- 
tions (Fig. 3b). Meanwhile, the electron density oscilla- 
tions observed in quasi ID case are absent in 2D case, 
even in the stationary state for the case with a bound- 
ary. In the latter case, we have also observed the onset 
of the luminescence ring at high enough pumping rate p 
(Fig. 4, cp. inset in Fig. 2a) and the shift of the ring posi- 
tion to larger r when increasing p. The result correlates 
well with the experimental dependence (Fig. la in Q) for 
pumping powers slightly above threshold. 

Qualitative picture. — The macroscopic electron-hole 
spatial separation and luminescence ring-shaped pattern 
formation can be explained as follows. At high photoex- 
citation power the carrier densities in the excitation spot 
quickly reach their critical values (~ d~ 2 ) when the re- 
pulsive in-layer Coulomb forces between particles become 
stronger than the attractive inter layer force. It leads to 
the appearance of in-layer electric fields ejecting the " sur- 
plus" electrons and holes from the excitation spot region. 
Then the e-h densities in the spot grow to the critical val- 
ues again and the ejection process recurs. Due to very 
high mobilities of the carriers in GaAs quantum wells 
at low temperature 18-2(J, the initial velocities of the 
ejected particles are also high. The absolute values of the 
velocities are restricted by the emission threshold of an 
optical phonon that gives u ma x ~ 10 7 cm/s [Hj]. in turn, 
the velocity ratio for the hot electrons and holes has the 
same order of magnitude as the ratio between their mo- 
bilities (fi e /fih > 10 for GaAs/AlGaAs QWs at tempera- 
tures T ~ 1 K [l9j.) Due to this, the electron distribution 
is mainly located much farther from the excitation spot 
than the hole one. During the in-plane expansion the 
electrons and holes cool down emitting acoustic phonons. 




FIG. 2: Ring-shaped luminescence patterns and spa- 
tiotemporal carrier dynamics in quasi ID case, (a) 

The in-layer distribution of luminescence for different critical 
velocities V c (main panel, p — 10) and for different pumping 
rates p in the excitation spot with size xq = 1 (inset, V c = 10). 
The distant peak of luminescence appears only at relatively 
small V c and the peak intensity grows with increasing p. Con- 
stant parameters are a = 0.2, c\ = 1, oz = 0.1, d 2 = 0.04, 
V e = 50, integration step At = 5 ■ 10 -4 . (b) Distributions 
of electrons (blue curve) and holes (red curve) for p — 10 
and V c = 10 (main panel). The luminescence distribution 
exhibits a peak at the overlap of electron and hole densities 
which are well separated. Electron density shows clear oscilla- 
tions (main panel and top-right inset). In the stationary state 
the creation rate of electrons and holes is balanced by the lu- 
minescence rate so the total number N of carriers saturates 
(bottom- left inset). 



In addition, their relative velocities decrease due to the 
intcrlayer Coulomb drag [22j. Finally, the luminescence 
ring appears in the overlapping region between the sta- 
tionary densities of the slowed electrons and holes. Note 
that at typical electron-acoustic phonon scattering time 
Te_ ac ~ 10~ 9 s [23] one comes to the promising estimate 
for the ring radius R ~ v max r e - a c ~ 0.1 mm (cp. [a|20j). 

Conclusion. — We have simulated a spatiotemporal 
evolution of the locally photoexcited electrons and holes 
localized in two separate layers. In contrast to previous 
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FIG. 4: The luminescence ring at different values of 
pumping rate p for 2D case with boundary (mirror 
reflection) rb = 100 in the stationary state. All other 
parameters are the same as for Fig. 3. 



observed drastically different electron and hole dynamics 
in ID and 2D cases and, more importantly, pronounced 
stationary spatial oscillations of the electron density in 
ID case. 

We thank A. S. Alexandrov, Yu. M. Kagan and F. V. 
Kusmartsev for helpful discussions. 



FIG. 3: Spatiotemporal dynamics of the luminescence 
ring and electron-hole distributions in 2D case. All 

spatial distributions are averaged over the polar angle, (a) 
Snapshots of luminescence distribution ni um (r, t) averaged 
over time window (ti/t s , t<z/t 3 ): (300, 400) for the red (I) 
curve, (400, 500) for the blue (II) curve, (500, 600) for the 
magenta (III) curve, (600, 700) for the brown (IV) curve. 
The peak of ni um mimics the ring-like patters observed in ex- 
periments 0, [H, HI, H3]. m some time after the pumping 
is switched on, the ring forms and expands even though the 
stationary regime is not reached (inset). Constant parame- 
ters: a — 1.5, excitation spot size ro = 4, p = 10, c\ — 1, 
c 2 = 0.25, d 2 = 2.25, V e = 50, V c = 10, At = 5 • 10 -4 . (b) 
The snapshots of electron and hole densities as well as lumi- 
nescence pattern averaged over time window (400, 500). The 
ring pattern forms at the overlap of legibly separated electron 
and hole distributions. 



theoretical studies we have shown that both 

in-plane spatial separation of electrons and holes and the 
ring-shaped luminescence pattern are formed essentially 
due to the in-layer Coulomb repulsion. A specific reason 
for the separation is found to be the difference in electron 
and hole mobilities. In turn, a specific reason for the lu- 
minescence ring is the dependence of exciton formation 
probability on the electron-hole relative velocity. We also 
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We have found a typo in the Y-projection of the interlayer electron-hole force in the MD code for 2D case. It 
affects the results shown in Fig. 3 and Fig. 4 of our original paper. Due to the typo the angular symmetry of carrier 
spatial distributions was broken that led to the absence of balance between pumping and decay in the 2D case without 
boundary and, eventually, to the permanent growth of the total number N of carriers in time t (inset in Fig. 3a). 

Simulations by the revised code have shown that there is a saturation of N(t) in the 2D case without boundary 
(see inset in Fig.l here). The balance occurs at about 350 particles that corresponds to the first plateau of N(t) 
shown on the inset in Fig. 3a. The spatial distributions of carriers and luminescence are quite similar to that shown 
in Fig. 3b, i.e., there is clear spatial separation of the electron and hole distributions though the luminescence ring is 
less pronounced (main panel in Fig.l here). 

The dependence of the ring position on the pumping rate p [l[ becomes very weak: the increase of p in few times 
leads to the growth of the ring radius on few percents Q. In this relation, the results shown in Fig. 4 are erroneous. 
Note that at p = 10 for the 2D case without boundary the maximal radius of the electron distribution is about 100 
so the boundary r& = 100 plays no role. 
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FIG. 1: Stationary spatial distributions of electrons (n e ), holes (nu) and luminescence (ni um ) averaged over time window (500, 
5000) in the 2D case without boundary. Inset: dependence N(t). All parameters are the same as for Fig. 3. 

So, the main consequence of the revised code is the saturation in the 2D case without boundary. Note that it brings 
an exciting opportunity to perform non-time-consuming molecular-dynamic simulations of spatiotemporal dynamics 
of charge carriers in mesoscopic 2D systems. 



[1] During MD simulation the time is changed by discrete steps, with the elementary time step At. The rate p was defined as 
a probability per unit time to create one e-h pair in the excitation spot, so that pAt < 1. One can refer this case as single 
generation regime. 

[2] It might result from that the pumping algorithm used [lj provides effectively too small pumping rate to apply the qualitative 
picture described in Sec. 5. 



